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ABSTRACT 

Simulation  of  the  nonlinear  effects  of  laser  atmospheric 
propagation  is  reported  in  this  paper.  A  set  of  atmospheric 
hydrodynamic  equations  and  the  nonlinear  electromagnetic  wave 
equation  for  interaction  with  the  atmosphere  are  given.  The  wave 
equation  describing  small-perturbation  instability  induced  by 
stimulated  thermal  Rayleigh  scattering  (STRS)is  derived.  Also, 
turbulence  is  theoretically  treated.  Based  on  these  equations, 
we  performed  a  set  of  4D  programs  (TURBLM)  to  simulate  thermal 
blooming  and  phase  compensation  for  it,  small-perturbation 
instability,  and  the  interaction  of  turbulence  with  thermal 
blooming.  The  simulated  results  are  in  agreement  with  theory  and 
experiment . 

Key  Words:  thermal  blooming,  small-perturbation  instability, 
turbulence . 


1 


1 .  Introduction 


The  propagation  of  a  high-power  laser  through  the  atmosphere 
is  an  extremely  complex  process,  even  without  considering  the 
nonlinear  process,  which  includes  a  sizable  number  of  factors. 
Among  these,  the  relatively  significant  nonlinear  effects  not 
involving  the  internal  processes  of  atmospheric  molecules  include 
thermal  blooming[l-3] ,  small -perturbation  instability  resulting 
from  stimulated  thermal  Rayleigh  scattering (STRS) [4-5] ,  as  well 
as  interaction  between  atmospheric  turbulence[6]  and  thermal 
blooming,  etc. 

Based  on  the  above-mentioned  scenario,  a  numerical 
simulation  was  conducted  in  this  paper,  and  a  4D  program  (TURBLM) 
was  written  to  investigate  the  foregoing  nonlinear  processes  of 
high-power  laser  atmospheric  propagation.  Apparently,  the 
experiment  on  laser  omnidistance  atmospheric  propagation  can  be 
well  simulated  if  a  linear  process  is  introduced  into  the  program 
without  the  necessity  of  making  significant  modification  to  its 
main  structure. 

2.  Thermal  Blooming 

During  high -power  laser  propagation  through  the  atmosphere, 
a  small  portion  of  energy  can  be  absorbed  by  the  molecules  and 
aerosol  in  the  atmosphere,  and  the  heated  air  can  lead  to  an 
increase  in  gas  pressure.  As  a  result,  the  gas  will  expand  at 
the  speed  of  sound  and  the  pressure  equalizes,  giving  rise  to  a 
decrease  in  density  and  lower  refractivity,  thereby  forming  a 
negative  lens  and  causing  laser  beam  divergence. 

When  there  is  a  laterally-blowing  wind,  the  air  density  of 
the  descending  wind  decreases,  bringing  about  lower  refractivity 
of  this  air  and  forming  a  unique  crescent  beam  bending  toward  the 
ascending  wind;  this  situation  can  produce  light  distortion. 
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bending  and  divergence,  eventually  affecting  light  beam  quality. 
This  kind  of  effect  is  called  thermal  blooming.  Thermal  blooming 
can  restrict  the  maximum  power  in  high-power  laser  propagation 
through  the  atmosphere,  and  this  is  one  of  the  most  critical 
problems  in  high-power  laser  atmospheric  propagation. 

There  are  two  basic  interacting  coupled  equations  in 
research  on  thermal  blooming,  namely,  the  light  wave  equation 
showing  the  modification  of  refractivity  caused  by  changes  in 
atmospheric  density,  and  the  hydrodynamic  equation  showing  the 
modification  of  atmospheric  density  change  resulting  from  laser 
heating. 


The  paraxial  wave  equation  is  as  follows: 

2il|S  +  VW  +  »'(^- 1)0-0 


(1) 


The  hydrodynamic  equation  under  constant  pressure 
approximation  is  as  follows: 

+  - 1^0/  (2) 

3t  a*  c; 

The  relation  between  atmospheric  refractivity  and  density  is 
as  follows: 

•  —  1  =  xp 

(3) 


From  Eq .  ( 3 ) , 

Se  —  ^  —  1  *=»  2(»  —  »#)/««  2kpi  ( 4 ) 

where  k  is  the  number  of  waves,  i|r  is  the  complex  amplitude  of 
optical  field,  i.e.,  E=4exp[i(kz-&t) ] ,  I  is  laser  light 
intensity,  i.e.,  I=(c/4ir)  |i|f  j^,  po  is  standard  atmospheric 
density,  Pi  is  laser  induced  perturbation  density,  p=po+pi,K  is 
the  Lorentz  coefficient,  v  is  wind  speed,  Cg  is  speed  of  sound, 
v=Cp/cVo.  Eqs.(l)  and  (2)  constitute  a  set  of  coupled  equations 
showing  the  interaction  between  laser  and  atmosphere,  which  form 
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the  foundation  for  research  on  the  linear  and  nonlinear  theories 
of  laser  atmospheric  propagation. 


Our  simulation  of  thermal  blooming  commences  from  this  set 
of  equations,  adopting  a  phase  screening  method  for  wave  equation 
(1)  as  well  as  the  fast  Fourier  transform  (FFT) . 


Suppose  iItr  is  a  field  at  2=Zr,  then  the  formal  solution  at 
Zj+j— Zg+Az  is 


•i+l 

f».+i  =  exp|  +  i'*  J  •  ^.| 


Under  the  binary  approximation,  the  following  equation  can  be 
formulated ; 

=  cxp((^Wi))exp(-|^<}i)cxp(^As*Vx)i». 

(6) 

where 

%+l 

and  let  8<l)i,=  (k/2)  and  let  it  represent  the  phase  distortion 

caused  by  thermal  blooming. 


It  can  then  be  seen  from  Eg.  (6)  that  owing  to  the  above 
treatment,  the  nonlinear  teim  is  simplified  and  acquires  a 
remarkable  physical  significance.  The  solution  to  Zg-Zg+i  can  be 
accomplished  in  three  steps;  (1)  Free  diffraction  of  the  optical 
field  Az/2;  (2)  The  nonlinear  phase  distortion  5(Pi,  introduced  due 
to  the  modification  of  atmospheric  refractivity  is  taken  into 
computation;  (3)  Repeated  free  diffraction  of  the  optical  field 
Az/2. 


Two  numerical  calculation  methods  are  adopted  for  Eg.  (2); 
one  is  FFT,  while  the  other  is  characteristic-type 
differentiation  method. 
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Based  on  Eqs.  (1)  and  (2),  we  edited  a  4D  program  (TURBLM) 
for  numerical  simulation.  The  result  of  the  simulation  is  in 
total  agreement  with  the  result  described  in  reference  8.  In 
addition,  another  simulation  was  also  made  for  thermal  blooming 
compensation,  which  resulted  in  a  great  improvement  in  light  beam 
quality[8] . 

3.  Small-perturbation  Instability 

Small -perturbation  instability  occurs,  during  stimulated 
thermal  Rayleigh  scattering ,  in  atmospheric  propagation  when  no 
adaptive  optics  is  used.  This  scattering  may  be  stimulated  by 
sound  waves  and  in  this  case,  it  is  called  stimulated  thermal 
Brillouin  scattering  (STBS).  Also,  it  may  be  stimulated  by  long¬ 
time  period  density  variation,  referred  to  as  stimulated  thermal 
Rayleigh  scattering.  The  relatively  long-time  period  STRS  was 
our  selection.  The  occurrence  of  small-perturbation  instability 
has  a  serious  practical  effect  on  high-power  laser  atmospheric 
propagation.  Atmospheric  turbulence  and  laser  noise  are 
stimulation  sources  for  small -perturbation  instability. 

Therefore,  for  high-quality  coherent  laser  beams  in  the  actual 
atmosphere,  one  must  suppress  the  instability  or  keep  the  growing 
small -perturbation  still  at  a  low  level. 

The  developmental  mechanism  of  STRS  instability  is  as 
follows:  laser  beams  with  small  perturbation  are  symbolized  by 
superimposition  of  unperturbed  intense  beams  and  weak  beams.  The 
intense  main  laser  (pvimp  light)  and  weak  perturbed  beam 
(scattering  light)  generate  interference  and  constitute 
interference  fringes.  The  interference  fringes  can  modulate  air 
absorption  and  heating  and  modify  the  refractivity.  The  intense 
light  can  enlarge  the  refractivity  area  and  transmits  energy 
further  to  the  weak  light  along  the  path,  thereby  intensifying 
it,  i.e.,  intensifying  the  perturbation. 
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There  are  two  methods  for  nvunerical  computation  of  the 
small-perturbation  instability.  One  is  to  subtract,  directly 
using  the  thermal  blooming  program,  the  optical  field  of 
perturbed  laser  atmospheric  propagation  from  the  optical  field  of 
unperturbed  laser  atmospheric  propagation.  The  result  of 
subtraction  is  the  temporal  and  spatial  variation  of  the 
perturbation  field.  This  method  is  simple  but  its  physical 
concept  is  not  clear.  The  other  method  is  as  follows:  when  the 
thermal  blooming  of  the  main  laser  (pump  light)  is  not  intense, 
it  is  possible  to  assume  that  the  main  laser,  not  subject  to 
atmospheric  effects,  undergoes  only  diffraction  and  absorption. 
And  the  modification  of  refractivity  is  basically  caused  by  the 
interference  from  the  main  laser  and  the  perturbation  field. 

The  perturbation  field  is  the  modification  and  increase  of  the 
scattering  light  in  the  atmosphere. 

Suppose  the  overall  optical  field  of  the  perturbation  field 
is  —  +  ,  and  the  air  density  variation  is  />;  =  p,  4- 

,  where  the  subscripts  p  and  s  are  quantities  related  to  pump 
and  scattering,  if'  and  p\  ,  respectively,  satisfy  Eqs.  (1)  and 
(2),  while  i|»p  and  Pp  also  satisfy  Eqs.  (1)  and  (2).  By 
subtracting  the  two  equations  of  i|r  and  p,  the  following  can  be 
derived : 


2ik  1^*  +  Vi^.  ==-  +  i>.p.,) 


(7) 


a  f,  +  t.ii!  ) 


(8) 


From  the  above  assumption,  p.,  I  ,  and  then 

the  simulation  coupled  equations  of  small-perturbation 
instability  are: 


Vl#,=  0 


(9) 
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(10) 


2it||‘  +  vl*. -"aVA 


3  p.  a  p. 


dt 


» —  1 


(11) 


Eq.  (9)  shows  the  free  diffraction  of  the  pump  light  in  the 
atmosphere  and  takes  only  atmospheric  absorption  into  account. 

Eq.  (10)  shows  an  increasing  wave  equation,  caused  by  the 
perturbed  optical  field  under  modulation  by  the  pumping  light  and 
density.  Eq.  (11)  shows  the  hydrodynamics  equation  for  the  case 
when  interference  by  the  piimp  light  and  the  perturbation  field 
causes  density  changes.  Eqs.  (9)-(ll)  constitute  a  set  of 
complete  coupled  equations  specifically  for  research  on  STRS 
induced  small -perturbation  instability. 


Among  the  simulation  parameters  that  we  adopted  are  the 
following;  the  light  intensity  of  the  initial  pvunp  field  is  a 
platform  distribution,  i.e.,  Ip=400W/cm^;  laser  wavelength  is 
10.6nm;  beam  diameter  D=50cm;  the  total  power  of  the  pump  field 
is  780kW;  the  light  intensity  of  initial  perturbation  is 
I.=0.04W/cm2;  the  period  of  the  perturbation  field  is  10cm; 
absorption  coefficient  a=0 . 2xl0-^cm-^•  and  horizontal  wind 
velocity  (along  the  x-direction)  v=10m/s. 


For  the  computed  perturbative  intensity  see  Fig.  1,  where 
the  curves  are  intensity  contours.  The  given  distributive 
fringes  of  initial  perturbative  intensity  form  a  45°  angle  with 
the  wind  direction  (along  the  x-axial  direction).  It  can  be  seen 
that  the  growth  direction  of  small  perturbative  fringes  (i.e., 
the  direction  of  the  contour  peaks)  gradually  approaches  the  wind 
direction  and  eventually  becomes  parallel  to  the  wind  direction. 
The  perturbative  field  first  breaks  up  and  gradually  turns  in  the 
wind  direction  before  connecting  one  with  another  to  form  fringes 
along  the  wind  direction.  As  a  result,  the  growth  direction  of 
the  perturbative  field  fringes  appears  parallel  to  the  wind 
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direction.  Such  kind  of  phenomenon  can  be  interpreted  like  this: 
since  a  small  segment  of  perturbative  fringe,  blown  by  the  wind, 
leads  to  the  decrease  of  air  density  in  the  descending  wind, 
while  the  light  is  diffracted  from  a  sparse  light  medium  to  a 
dense  light  medium,  this  small  segment  of  perturbative  fringe 
approaches  to  the  wind  direction.  When  the  propagation  distance 
and  time  are  sufficiently  long,  the  perturbative  fringes  are 
parallel  to  the  wind  direction. 

Additionally,  with  the  growth  and  change  of  perturbation  in 
space,  its  perturbative  period  generally  does  not  change.  In 
Fig.  1,  the  number  of  fringes  stays  at  5  constantly.  But  it 
seems  that  there  are  10  fringes  at  the  transmitting  aperture. 

This  is  because  the  sinusoidal  distribution  is  both  negative  and 
positive.  Actually,  there  are  only  5  fringes.  The  above  two 
arguments  conform  to  the  SABLE  experiment  (Scaled  Atmospheric 
Blooming  Experiment )  as  described  in  reference  9 .  The 
calculations  in  the  latter  are  listed  in  Fig.  2,  where  the 
simulation  parameters  are:  number  of  thermal  blooming  distortion 
Nd=220;  wind  fluctuation  a=ar/v=0.45  (a  is  variance);  niimber  of 
perturbation  modulation  is  5;  transmitting  aperture  is  35cm,; 
wavelength  is  2.9^m;  power  is  lOkW;  transmission  distance  z=400m. 

Fig.  3  shows  the  development  of  perturbative  intensity  with 
transmission  distance  z  at  two  points  (a:  x=y=0;  b:  x=20cm,  y=0) 
in  horizontal  space.  The  growth  and  variation  of  small 
perturbation  can  be  seen  in  this  figure:  the  perturbation  tends 
to  grow  in  space  when  t=0.005s.  It  grows  extremely  rapidly 
because  damping  is  not  taken  into  consideration.  For  instance, 
the  spreading  term  of  refractivity  and  random  wind  are  not 
considered  to  suppress  the  growth  of  small  perturbation. 
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x(cm)  x(cfn) 


Fig.  1.  Perturbative  intensity  contour  along  z-direction 
at  t=0.005s 

(a)  z=0  (b)  z=lkin  (c)  z=2kin  (d)  =3kin 


direction  of  wind 


Fig.  2.  Intensity  contour  on  the  beacon 


9 


Fig.  3.  Perturbative  intensity  with  propagation 

distance  Z  (a)  J  _  ^  a.  0|  (t)  x  •  20  cm,  jr  »« 0 

Let  us  briefly  analyze  Fig.  3.  According  to  Karr's  thermal 
blooming  linear  analysis[5],  the  amplitude  of  perturbative  field 
intensity  changes  in  time  and  space  variation  as  follows: 

/.(=,<)  ocexp[3.T(^)''*)  (12) 


where  NA=rk2t/2Tt  is  the  number  of  thermal  blooming  waves; 
r=(8.25xlO‘^cmVJ)aIp;  =  is-  the  perturbation  Fresnel 

number;  kp  is  the  number  of  perturbative  waves.  When  t=0.005s, 
by  substituting  the  simulation  parameters  in  Eg.  (12), 

/,(c,0  oc  expC3.  { 13  ) 

where  z  is  in  km.  From  the  imaginary  line  in  Fig.  3,  it  can  be 
calculated  that  (a)  Ig  varies  as  exp[ 4 .  ,  (b)  I3  varies  as 

exp[ 2 . 8z^^^] .  Therefore,  we  know  that  the  calculations  do  not 
closely  agree  with  the  analysis  results  in  Eg.  (13).  This  is 
because  Karr's  analysis  did  not  consider  the  effect  in  the 
horizontal  direction  but  the  horizontal  infinite  plane  wave, 
while  in  our  calculations,  the  effect  in  horizontal  direction  is 
taken  into  account.  As  a  result,  the  growth  of  perturbative 
intensity  in  horizontal  direction  is  different  in  both  cases. 
However,  our  calculations  generally  conform  to  linear  eguation 
(13).  Besides,  when  the  perturbative  field  grows  to  a  certain 
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degree,  its  growth  rate  will  slow  down  as  shown  in  Fig.  3(a). 
4.  Effect  of  Atmospheric  Turbulence  on  Thermal  Blooming 


The  variation  process  of  atmospheric  turbulence  is  a  random 
process,  which  is  to  be  described  with  a  statistical  law.  We  use 
the  Karman  spectrum  as  its  statistical  law.  Atmospheric 
turbulence  has  two  kinds  of  effects  on  theimial  blooming,  that  is; 
the  turbulent  beam  expansion  can  cause  changes  to  intensity 
distribution  and  thus  directly  affect  thermal  blooming,  and 
turbulence  perturbation  can  affect  temperature  gradient  and  cause 
thermal  blooming  to  change.  The  computation  of  the  combining 
action  of  turbulence  and  thermal  blooming  is  similar  to  the 
computation  method  without  turbulence,  but  unlike  the  latter,  it 
is  required  to  add  the  phase  shift  6(pt  introduced  by  turbulence 
and  the  phase  shift  Sq)b  produced  by  thermal  blooming  to  form  a 
phase  shift  caused  by  interaction  between  turbulence  and  thermal 
blooming  S(p=6(Pt+6(Pb.  The  phase  distortion  created  by  turbulence 
is 

*•+1 

<5^  =  ^  J* 


where  the  subscript  t  represents  turbulence  as  a  cause.  Thus, 
the  following  equation  can  be  derived  to  calculate  de,(*,y)  . 


-f-OO  4’®0 


de,Cx,y)  =  J  J  di.di',«xpCt(i.*  + 


(14) 


where  is  the  Karman  spectrum 


<P.(10  =  0.  033Cl(iS  +  ij  +  if  +  if)  ( 15  ) 


where  Cf  is  the  refractivity  structure  constant,  ko=2n:/Lo,  where 
Lo  is  the  external  scale  of  turbulence.  In  this  paper,  Lo=10m; 
kx,  ky,  and  kj,  are  components  of  wave  vector;  and  a(k2,ky)  is  the 
two-dimensional  complex  Gaussian  random  function. 


11 


By  making  the  integration  of  Eg.  (14)  discrete,  the 
following  can  be  derived: 


1^9, z)  — 


•  Ca'(w»,i«)  +  to*(ni,«)l 


where  a' (m,n)  and  a"(m,n)  are  the  Gaussian  distribution  functions 
with  variance  1;  N  is  the  number  of  sample  points;  Ax  and  Ay  are 
the  intervals  between  sample  points;  kAx  and  lAy  are  the  space 
location  coordinates.  Eq.  (16)  is  a  working  formula  for  phase 
distortion  caused  by  atmospheric  turbulence. 


Using  Eq.  (16),  we  computed  the  interaction  between 
turbulence  and  thermal  blooming  with  the  following  simulation 
parameters:  the  incident  beam  is  a  Gaussian  beam;  laser 
wavelength  p=10.6iim;  power  p=80kW;  focal  length  f=6km;  beam 
diameter  D=50cm;  absorption  coefficient  a=0 . 2xl0‘^cm'^;  horizontal 
wind  velocity  v=10m/s. 


Fig.  4  (a),  (b),  (c)  and  (d)  respectively  are  intensity 
contours  on  the  focal  plane  when  Cl-0,  10*®m"^^^,  5xl0*^*m'^'^, 

where  cf-0  is  a  pure  thermal  image  without  turbulence. 

It  can  be  seen  that  turbulence  has  a  rather  great  effect  on 
thermal  blooming.  The  more  intensive  the  turbulence  (the  greater 
the  Oi),  the  greater  the  laser  beam  expansion  and  trembling.  It 
can  even  cause  beam  to  expand,  resulting  in  further  decrease  of 
beam  quality. 
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Fig.  4.  Intensity  contour  at  focal  plane 
with  thermal  bloomino  and  turbulence 

(«)  Cl  -  0|  ft)  Cl  - 

(c)  Cl  •  5  X  10-‘*  (rf)  Cl  —  10-'»  m-*^ 

When  there  is  no  turbulence,  the  number  of  the  thermal 
blooming  distortion  is 


/oas^dw/dr  ( 17  ) 

RoPoc^a 


Based  on  numerical  simulation  and  experimentation,  Gebhardt 
obtained  the  following  formula: 

, _ 1 

-  1  +  0.  0625Ar*  (18) 

where  Irei  is  the  ratio  between  the  intensity  with  thermal 
blooming  and  intensity  under  linear  absorption,  which  is  a 
parameter  indicating  the  beam  quality;  when  Ne=0,  Irei=l/  i.e., 
the  beam  is  not  distorted.  By  substituting  simulation  parameters 
in  Eg.  (17),  No=23,  and  Irei=0.03  can  be  obtained.  Based  on  Fig. 
4(a),  Irei=0.035  can  be  computed,  which  is  close  to  the  calculated 
value  from  Eg.  (18).  For  images  with  turbulence,  1^^1=0.028, 
0.015,  and  0.012  can  be  calculated,  respectively,  from  Fig.  4(b), 
(c),  and  (d).  It  is  known  therefore  that  turbulence  has  a  small 
effect  on  thermal  blooming  when  C*  =10'^*m’^'^,  but  it  has  a  strong 
effect  on  thermal  blooming  when  5x1 
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5.  Conclusions  and  Discussion 


The  following  conclusions  can  be  reached  based  on  the  above 
analysis  and  calculations  iri  this  paper: 

( 1 )  There  is  small-perturbation  instability  in  laser 
atmospheric  propagation,  which  has  an  unfavorable  effect  on  the 
main  laser; 

(2)  Small-scale  growing  fringes  approach  toward  the  wind 
direction,  and  the  perturbation  period  basically  remains 
unchanged ; 

(3)  The  strength  of  atmospheric  turbulence  has  different 
effects  on  laser  atmospheric  propagation. 

This  paper  mainly  discussed  thermal  blooming,  STRS  induced 
instability  and  the  interaction  between  theinnal  blooming  and 
turbulence.  It  is  noteworthy  that  there  is  still  small- 
perturbation  instability  caused  by  the  stimulated  thermal 
Brillouin  scattering  (STBS)  in  high-power  laser  atmospheric 
propagation.  The  major  difference  between  STBS  and  STRS  lies  in 
that  there  is  frequency  shift  between  scattering  light  and 
incident  light  in  STBS  {k*0 ,  Since  STBS  takes  place  in  the 

transient  process  of  thermal  blooming  and  the  transient  process 
is  only  3-5  times  the  hydrodynamics  time  tH(tH«D/C3;  D  is  beam 
diameter,  is  the  speed  of  sound),  STRS  has  a  small  effect  on 
the  continuous  laser  beams.  However,  its  effect  on  pulsed  laser 
propagation  cannot  be  neglected. 

This  paper  was  received  for  editing  on  March  4,  1993.  The 
edited  paper  was  received  on  June  27,  1995. 
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